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Abstract 

In a previous paper [Phys. Rev. 50, 1980 (1994)] we gave account of the 
nondispersive band first found experimentally at -4.4 eV for CdTe(OOl) by 
Niles and Hochst. We have characterized this band as a surface-induced bulk 
state. In a second paper we showed that a similar state does exist in II- 



. VI and III-V zincblende semiconductor compounds. In this paper we show 

I . that there are more such states within the valence band energy interval. We 

^ . use tight-binding hamiltonians and the surface Green's function matching 



method to calculate the surface and surface-induced bulk states in the wide 



^ . band gap zincblende semiconductors CdTe, CdSe, ZnTe and ZnSe. We find a 

distinctive surface state for the cation and two for the anion termination of 
the (OOl)-surface and three (OOl)-surface-induced bulk states with energies 
that correspond to the value of the heavy hole, light hole and spin-orbit bands 
at X. 

PACS:73.20.At; 73.20.Dx 
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I. INTRODUCTION 



The study of the physics of surfaces, interfaces, quantum wells and superlattices of 
semiconductors has call interest in the last few years. [jlj-Q The interest is not only for 
binary compounds but also, more recently, for ternary and quaternary compounds. 

The starting point is an accurate description of the electronic band structure of the binary 
compound by a method which serves as a sound basis to a clear and simple description of 
the more complicated compounds and systems. 

In previous work |TD|,|TT| we have used the tight-binding formulation to calculate the 
total density of states (DOS) and, in conjunction with the known surface Green's function 
matching method (SGFM), [T^ to derive the surface and the interface DOS. In this paper we 
want to use the same method to calculate the surface and surface-induced bulk states. The 
method can be also applied to superlattices P,pr3| and phonons. These band structures 
can be also used to calculate transport properties in heterostructures as quantum wells, for 
example, by making use of the well known many-body formulation by Keldysh. [|1^] 

In this paper, we consider the band structure of the II- VI wide band gap semiconductors 
CdTe, ZnTe, CdSe and ZnSe. We have first obtained the bulk bands (infinite medium) 
from the direct diagonalization of the tight-binding Hamiltonians and from the poles of 
the real part of the corresponding Green's function which is totaly equivalent. These very 
same results were obtained using the SGFM method, from the (OOl)-bulk-projected Green's 
function. This is actually a proof of consistency which gives us confidence in the new results. 

Our calculation shows the existence of three surface-induced bulk states which have no 
dispersion from F to X and which do not change in energy with different atom termina- 
tion of the sample. These states were first found experimentally by Niles and Hochst in 
CdTe(OOl). There is also a distinct surface state for each the anion and the cation termi- 
nation of the surface. They do differ noticely in energy. A general pattern develops for the 
location of the different bands in II- VI wide band gap zincblende semiconductors within the 
range of their valence band energy. 
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The rest of the paper is organized as follows. In Section II, to make the paper self- 
contented, we briefly describe the main highlights of the method and present the formulae 
that we have used. In Section III, we first introduce the general characteristics of the valence 
band and the location in energy of the surface and surface-induced bulk states. Next, we 
discuss the results for each system studied. In a final section, we present our conclusions. 
We have included an appendix where we quote all the tight-binding parameters of interest 
for this calculation. 



II. METHOD 

We make use of tight-binding Hamiltonians. Since the Green's function matching method 
takes into account the perturbation caused by the surface exactly, at least in principle, we 
can use the tight-binding parameters (TBP) for the bulk. |10|,|ll|,|16| This does not mean 
that we are using the same TBP for the surface and the bulk. Their difference is taken into 
account through the matching of the Green's functions. We use the method in the form 
cast by Garcia-Moliner and Velasco. |T^ They make use of the transfer matrix approach 



first intoduced by Falicov and Yndurain. |]T^ This approach became very useful due to the 



quickly converging algorithms of Lopez-Sancho et al. ||T8| Following the suggestions of these 
authors, the algorithms for all transfer matrices needed to deal with these systems can be 



found in a straightforward way. ||T9| This method has been employed successfully for the 



description of surfaces, P,|T0|,|T6| interfaces, ||TT|,|20[| and superlattices. 0,1 



A. The formalism 



We have first, calculated the bulk (infinite medium) band structure of the compounds by 



the tight-binding method (TB) in the Slater-Koster language using an orthogonal basis 
of five orbitals, sp^s* . The s* state is introduced to properly locate in energy the conduction 
band usually formed by d states in the II- VI zincblende (ZB) semiconductor compounds. 
2^ , |23[| In our calculation, we have included the effect of the spin-orbit (SO) interaction. ||2^ 



The TBP that we have used in our calculation are listed in the appendix. They reproduce 
the known bulk bands quite well. P J25|J26|] We assumed ideal truncation. 



We obtain the Green's function from 

{uj-H)G = I (2.1) 

where u is the energy eigenvalue, H the tight-binding Hamiltonian and / is the unit matrix. 
We adopt the customary decription in terms of principal layers. We label them with positive 
numbers and zero for the surface principal layer. Atomic layers are labelled with negative 
numbers and zero for the atomic surface layer. Let | n) be the principal wave function 
describing the rfi^ principal layer. It is a LCAO wave function with one s-like, three p— like, 
and one s*-like atomic functions per spin on each atom in the unit cell (there are two different 
atoms in the cell, and two atomic layers per principal layer, i. e., it is a 20-dimensional 
vector). If we take matrix elements of eq. ( |2.1[ ) in the Hilbert space generated by the 
complete set of the wave functions | n), we get 

{n\{uj-H)\m) = 5mn. (2.2) 

Since there is only nearest-neighbors interactions between principal layers, the identity 
opertator for the n^^ principal layer is 

/ =1 n - 1) (n - 1 I + I n) (n I + I n + 1) (n + 1 I (2.3) 

and therefore Hm,m+i = for i > 2. By inserting eq. ( |2.3|) into eq. (|2.2|) we get 

The matrix elements of the Hamiltonian, Hnm-, that appear in this formula are 2x2 super- 
matrices (each principal layer contain two atomic layers) each of whose elements is a 10x10 
matrix (since we are using a 5-wave functions basis by spin per atom). For example, for the 
surface 



00 



^ /loo ho^i ^ 



(2.5a) 



01 



V 



(2.5b) 



ho-2 /iO-3 

h-i-2 

Notice that rows are labeled with the index of the surface principal layer zero (containing 
atomic layers and -1) while the columns are indexed with the first principal and second 
principal layer (atomic layers and -1 and atomic layers -2 and -3), respectively. We shall 
adopt the hypothesis of an ideal, non reconstructed surface. Then, for the (OOl)-surface, for 
example, we have one atomic layer of anions and one of cations per principal layer. Therefore 
in this case /iqo 7^ ^-i-i, ^o-i = ^-lo- Thus, to calculate Hqq and Hqi we need to know 
^00) ^-1-1, /lo-i, and /i_i_2- These matrices are readily written in a tight-binding language 
and can be calculate with the bulk parameters as mentioned above. They depend on the 
wave vector k. 

Using eq. ( |2.4| ) for n = m, and m = for the surface, it is straightforward to get the 



surface Green's function [12 



G:' = (uI-Hoo)-Ho,T 



(2.6) 



and the principal layer projected bulk Green's function [O 



It is customary to define the transfer matrices as 



(2.7) 



G 



k+lp 



TGkp, G 



k+lp ^kp 



GkpS, k > p >0 



Gjj-fi TGij, 



G 



ij+i 



GijS, j>i>0. 



(2.8a) 



(2.8b) 



These matrices can be calculated by the quick algorithm of Lopez-Sancho et al. |jT8[, recal- 
culated later by Baquero |T9| (see Refs. p^p!T| for a compilation of all the formulae). 

From the knowledge of the Green's function, the surface states and the surface-induced 
bulk states can be calculated from the poles of the real part of the corresponding Green's 
function. We have applied previously this formalism to other surfaces [^,[TD|,|TB|,|2^ . We will 
now present our results. 



III. RESULTS 



The general characteristics of the zincblende II-VI wide band gap semiconductor valence 
band as we have obtained from our calculations is as follows. The heavy-hole (hh) and light- 
hole (Ih) bands follow each other closely in energy. The corresponding wave functions are 
mainly [p^, Py) in character. The hh band disperses from F to X about 2.0 eV and the Ih 
one about 2.4 eV. The spin-orbit splitting is around 1 eV in the Te-compounds and around 
0.5 in the Se-ones. This band reaches X at about 5.0 eV. This is the band with the most 
dispersion and is composed essentially of Pz states only. Finally, a very deep bulk band of 
mostly s— character appears below -10 eV. It disperses about 2 eV from V — X. We will call 
this band biQ. 

From the (OOl)-bulk-projected Green's function we get also the energy of the (001)- 
surface-induced bulk states. Three such states appear, B^, Bi and B^. These surface- 



induced bulk bands show no dispersion |5,2q,2£]. B^ was first found experimentally by Niles 



and Hochst and confirmed later by Gawlik et al. |^ for GdTe(OOl). Bh mixes with the hh 
band at X and is located at the same energy (in F) than the surface state, (S'c, see below) of 
the cation terminated (OOl)-surface. Bi mixes with the Ih band at X. The composition of 
both Bh and Bi is mainly of (jo^., Py), while the one of B^ which mixes with the spin-orbit 
band at X is (s, Pz)- The three states appear at the same position in energy irrespective of 
the cation or anion termination of the surface as one expects for surface-induced bulk states 
which only depend on the surface through the boundary condition (the wave function has 
to be zero at the surface). 



The (OOl)-surface valence band is rich [g,^] in other features. In particular, three 
characteristic surface states do exist in this range of energy. Two correspond to the anion 
{Sal, Sa2) and the other to the cation (S'c) termination of the (OOl)-surface. In all the 
systems considered, the anion terminated surface higher band {Sai) follows roughly the 
dispersion of the hh bulk band but it is at slightly higher energy. The cation terminated 
surface band (S'c) starts roughly around 2-3 eV from the top of the valence band in F and 
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has a varying ammount of dispersion. The two states appear at very different energy values 
and are distinctive of the termination of the surface for the four systems under consideration. 
One could speculate on the use of these two states to characterize the termination of the 
surface. A second (OOl)-surface anion state {Sa2) appears at much lower energies, for all 
the compounds considered, near the lowest bulk band 610. This surface state has very little 
dispersion and its wave function is mostly s— character. 

There is only one set of experimental results available at this moment to compare our 
results with, namely the already cited CdTe(OOl). We have studied this state and 
found that this is a surface-induced bulk one. Such bulk states appears as a consequence 
of the breaking of the symmetry implied by the creation of the (OOl)-surface in the infinite 
medium. To support our argument, we showed IQ specifically that a pole in the real part 
of the (OOl)-bulk-projected Green's function appears at -4.4 eV from F to X, and does not 
appear neither in the real part of the (OOl)-surface-projected Green's function nor in the 
real part of the bulk Green's function for any termination (anion or cation) of the sample. In 
other words, this is neither a solution of the Shrodinger's equation for an infinite medium nor 
a solution for a semi-infinite one that is localized at the surface. It is a bulk state existing as 
a consequence of the creation of the surface for which the crystal momentum is not a good 
quantum number anymore due to the new boundary condition. It has no e(A;)-dispersion as 
a consequence. This is a surface-induced bulk state, as we already stated above. 

In the rest of this paper, we want to describe in detail these (OOl)-surface and (001)- 
surface-induced bulk states for the four systems under consideration as well as the bulk 
bands obtained from our calculation. We will use our Hamiltonians to describe interfaces, 
superlattices and quantum wells of these compounds in future work. 

We summarize our results in Tables I-IV. In Table I, we give the characteristics of the 
bulk bands according to our obtained values. We give the gap and the spin-orbit coupling 
for each element and the value for the bands in F and at X. In Table II, we give the energy 
of the three surface-induced bulk states, then, in Table III, we give the energy of the surface 
states and, finally, in Table IV, the composition of all the bands involved in our work. 
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A. CdTe(OOl) 



Fig.l shows the full valence band for this system. The bulk states appear as the dotted 
curves in the figure. Our tight-binding description reproduces correctly the band-gap value 
(1.6 eV) and the spin-orbit splitting (0.9 eV). Our tight-binding parameters were adjusted 
to give energy values at the X-high symmetry point of the Brillouin zone that reproduce 
the experimental values for CdTe very closely. The heavy hole and light hole bands 

have a width of 1.8 eV and 2.2 for this interval of the Brillouin zone, respectively. The wave 
functions are constituted mainly by an admixture of Px and Py states. See Tables I and IV. 

The surface states existing in this energy interval are identified as Sai, Sa2 and Sc- The 
band Sai corresponds to the anion termination of the surface (black triangles), its width 
(r — X) is about 1 eV. The wave functions are constituted mainly by an admixture of s 
and Pz cation orbitals. The Sc surface state (white triangles) corresponds to the cation 
terminated (OOl)-surface. This band develops from about 2.2 eV below the top of the 
valence band in F to about -3.5 eV at X. Its width is therefore about 1.3 eV and is mainly 
of s-character. This band mixtures with a surface-induced bulk band, {Bi)^ for about half of 
the wave- vector interval T — X (see below). The Sa2 surface state shows a small dispersion 
and is located at -8.5 eV. Its wave function is of s— character. These surface states were 
described before. |^ 

Three surface-induced bulk states exist in this range of energy according to our findings. 
They are denoted by Bh^ Bi and Bg. 

The states B^ and Bi are newly found surface-induced bulk states. Bh is located at 
-1.8 eV and has no dispersion as is characteristic of these surface-induced bulk bands. B^ 
seems to mix with the hh band at X. Bi appears at -2.2 eV and mixes with the Ih band 
at X. Bi could present branches with slighly different energy and be actually a group of 
states differing slightly in energy. This states are mainly {p^, Py). Bi starts at F at the 
same energy as the surface state Sc- Their energy difference is very small for about half the 
interval T — X when the surface state begins to have an important dispersion while the bulk 
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one has none. We have obtained each state from the corresponding Green's function. The 
Bg state has a (s, p2)-composition. 

In conclusion, CdTe(OOl) presents in the range of energy covered by its valence band, 
as it is shown in Fig. 1, in adition to the known bulk bands distinct surface states for 
each the anion {Sai, Sa2) and the cation (Sc) terminated (OOl)-surfaces and a total of three 
surface-induced bulk states, Bh, Bi, and Bs- Bs is the recently stablished surface-induced 
bulk state at -4.4 eV. None of them presents dispersion and they appear at the same 
energy irrespective of the termination of the (OOl)-surface. 



B. CdSe(OOl) 



In general, CdSe is grown in an hexagonal structure. Nevertheless, recent experiments 



3T| have shown that it is possible to grow this material in a cubic phase and that it is stable. 
Our calculation has been done assuming the zincblende structure. Our bands do not apply 
for the hexagonal structure. 

We can see in Fig. 2 the general characteristic of the II-VI zincblende wide band gap 
semiconductor valence band already described above. The hh and Ih bands follow each other 
closely in energy and are both composed by {px, Py). The hh band disperses from F to X 
2.2 eV and the Ih one 2.4 eV. The spin-orbit splitting is 0.43 eV and the band (pz) ends at 
-5.0 eV at X and has therefore a width of about 4.6 eV. The feio band appears at -11.0 eV 
in F and has little dispersion towards X. 

For the anion terminated (OOl)-surface, there is one surface state (Sai) of (s, pz)- 
composition which develops right from the top of the valence band in F to X with a quite 
less dispersion than the hh band. At X it takes the value 1.0 eV. The cation terminated 
(OOl)-surface shows a surface state (Sc) of s-character at a lower enery in F (-2.2 eV) that 
disperses very little up to half the T — X interval but then disperses strongly and ends at 
-4.0 eV at X. Its width in the interval is 1.8 eV. The Sa2 band appears again very closely to 
the bio bulk one and shows very little dispersion as it is the general behavior for this surface 
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band. 

The three (OOl)-surface-induced bulk states {B^, Bi and Bg) appear at -2.2 eV {px, Py), 
-2.45 eV {px, Py) and -5.0 eV (s, pz) irrespective of the cation or anion termination of the 
surface. The first two mix with the hh and Ih bands at X, respectively, and the third one 
mixes with the spin-orbit band at X. None of them presents dispersion. 

As we can conclude from the previous description, CdSe(OOl) follows very closely the 
general picture described above for the number, composition and relative location of the 
surface and surface-induced bulk states within the valence band energy interval. We will 
now show that the last two compounds studied follow the same pattern. 



C. ZnTe(OOl) 

The full valence band for ZnTe(OOl) is presented in Fig. 3. It is clear that the general 
pattern is followed. Sai is the highest band in energy. It starts at F and disperses 1.0 eV as it 
reaches X. It is of (s, Pz) character. The hh and Ih bands (both [p^, Py)) follow each other 
and disperse 2.0 and 2.3 eV, respectively. B^ and Bi both mix with the hh and Ih band, 
respectively, at X and have the same composition. The Sc state disperses 2.0 eV and mixes 
with the Bh band at F. It is s-like character following the general rule. Finally, the Bs state 
is located at -5.3 eV from the top of the valence band and is of (s, ]92)^character. This state 
seems to be composed of two very nearly lying in energy states. The surface-induced bulk 
states present no dispersion and are found at the same energy irrespective of the surface 



termination. |28| There is no experimental evidence known to us for this findings. The two 



surface states were described before. 



D. ZnSe(OOl) 

We present our results for this compound in Fig. 4. In all respects ZnSe(OOl) follow 
the general picture for the states in the valence band energy interval. The hh, Ih and 
spin-orbit bands have the usual composition and behavior. Three different surface states 
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appear, two for the anion-surface (S'ai, Sa2) and one for the cation-surface (Sc), and three 
surface-induced bulk states with the same usual characteristics. The surface states have 
been studied previously. |^ We have summarize the details in Tables I-IV. 

IV. CONCLUSIONS 

In conclusion, we have studied the (OOl)-surface and surface-induced bulk states that 
appear in the range of the valence band energy for the II-VI wide band gap zincblende 
semiconductors. We found that a general pattern applies to the states and relative energy 
of these states with respect to the bulk bands. First, near the top of the valence band a 
surface state develops with little dispersion from F to X for the anion-terminated surface. 
A surface state also exists for the cation terminated surface but it appears at a quite lower 
energy. A second anion-surface state appears in the range of energies of the botton of the 
valence band. Simultaneously three (OOl)-surface induced bulk states appear which show 
no dispersion and appear at the same energy irrespective of the surface termination. These 
states mix with the hh, Ih and spin-orbit bands at X, respectively. 

It would be interesting to study other directions of the surface as well as interfaces, 
superlattices and quantum wells in different high-symmetry directions to characterize the 
different states that this break of symmetry due to the creation of a border might induce. 
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APPENDIX: 

We present in this appendix the tigh-binding parameters that we used for the systems 
studied in this work. 

Tight-binding parameters used in our calculation. We used the notation of Bertho et al. 



Parameter 


CdTe^ 


ZnTe'^ 




ZnSe^ 


CdSe'^ 




-8 


19210 


-9.19000 


-12 


42728 


-10.16740 


e; 





32790 


0.62682 


1 


78236 


1.03400 


E's 


-0 


95000 


-1.42000 





04728 


1.07977 


Ep 


6 


93790 


3.77952 


5 


52031 


7.64650 


Vss 


-5 


00000 


-6.64227 


-6 


50203 


-2.89240 




2 


iobOO 


1.94039 


3 


30861 


3.01320 


Vxy 


4 


52817 




c; 
o 




^ 7'?nzLn 

o. t ou^u 


yac 
' sp 


3 


31200 


5.92472 


1 


13681 


2.16040 


yea 
' sp 


3 


63824 


4.67265 


5 


80232 


5.65560 


VI * 


10 


44540 


6.22682 


7 


84986 


6.02650 


V", * 

' S*S* 


6 


62960 


6.77952 


8 


52011 


3.96150 


VT 

s*p 


2 


52468 


2.96202 


3 


26633 


2.11640 


Vf 

y gtp 


2 


94540 


3.82679 


1 


86997 


2.21680 


Xa 





32267 


0.36226 





19373 


0.14300 


Ac 





07567 


0.02717 





01937 


0.06700 



^ Ref. § 
^ Ref. H 
^ Ref. M 
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FIGURES 

FIG. 1. The full electronic valence band structure for CdTe(OOl) in the T — X direction. The 
dot lines are the bulk band structure from direct diagonalization of the tigh-binding Hamiltonian 
using the tight-binding parameters that appear in the Appendix. The two surface states for the 
anion terminated (OOl)-surface are denoted by full triangles, Sai and Sa2- The empty triangles 
are the surface state for the cation terminated surface, Sc- The (OOl)-surface-induced bulk states, 
Bfi, Bi and Eg are the full dots (the dashed lines are keep up only for eye guide). Their occurence 
is explained in the text. Bg was first found experimentally by Niles and Hochst.[3] 

FIG. 2. Electronic structure of the CdSe(OOl) valence band in the T — X direction. The 

conventions are the same as in Fig. 1 (Sec figure caption). The (OOl)-surface-induced bulk states, 
Bfi, Bi and Bg mix at X with the hh, Ih and spin-orbit bands respectively as a general rule. The 
pattern presented by the surface and surface-induced bulk states as they appear in this figure is 
the general one found also in the rest of the systems studied. 

FIG. 3. Electronic structure of the ZnTe(OOl) valence band in the V — X direction. See figure 
captions 1 and 2 for details. 

FIG. 4. Electronic structure of the ZnSe(OOl) valence band in the V — X direction. See figure 
captions 1 and 2 for details. 
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TABLES 

TABLE I. Characteristics of the band structure. Eg is the gap and Ego is the spin-orbit 

splitting. The last three columns arc the values at the X point of the Brillouin zone of the energy 
of the heavy hole band (Ehh), the hgh hole band (Eih), and the spin orbit {Ego)- All the values 
are in eV. 





X-point 






Eso 


Ehh 


Eih 


Eso 


CdTe 


1.602 


-0.9 


-1.7 


-2.2 


-4.4 


CdSe 


1.78 


-0.4 


-2.18 


-2.36 


-4.89 


ZnTe 


2.39 


-0.91 


-1.93 


-2.40 


-5.50 


ZnSe 


2.82 


-0.45 


-1.95 


-2.19 


-5.30 



TABLE II. The energy position of the (OOl)-surface-induced bulk states. All the values are 

in eV. 

CdTe CdSe ZnTe ZnSe 

Bh -1.8 -2.2 -2.0 -2.0 

Bi -2.2 -2.45 -2.5 -2.3 

Bs -4.4 -5.0 -5.5 -5.3 
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TABLE III. The energy position of the (OOl)-surface states. All the values are in eV. 



F-point 


X-point 




Sal 




So 


Sal 


Sa2 


So 


CdTe 


0.0 


-8.5 


-2.2 


-1.0 


-8.5 


-3.6 


CdSe 


0.0 


-10.4 


-2.2 


-0.9 


-10.4 


-3.9 


ZnTe 


0.0 


-11.6 


-2.0 


-0.9 


-11.7 


-4.4 


ZnSe 


0.0 


-12.5 


-2.0 


-0.7 


-12.4 


-4.2 



TABLE IV. The wave function decomposition for all the states that appear in this study. 



State Composition 

hh {Px,Py) 

{Px,Py) 

so {pz) 

Sal {s,Pz) 

Sa2 is) 

Sc (s) 

Bh {Px,Py) 

Bl iPx,Py) 

Bs {s,Pz) 
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